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ABSTRACT 

Context. The determination of ages of central stars of planetary nebulae (CSPN) is a complex problem, and there is 

presently no single method that can be generally applied. We have developed several methods to estimate the ages of 

CSPN, based both on the observed nebular properties and in some properties of the stars themselves. 

Aims. Our aim is to estimate the ages and the age distribution of CSPN and to compare the derived results with mass 

and age determinations of CSPN and white dwarfs based on empirical determinations of these quantities. 

Methods. We discuss several methods to derive the age distribution of CSPN, namely, (i) the use of an age-metallicity 

relation that also depends on the galactocentric distance, (ii) the use of an age-metallicity relation obtained for the 

galactic disk, and (iii) the determination of ages from the central star masses obtained from the observed nitrogen 

abundances. 

Results. We consider a sample of planetary nebulae in the galactic disk, most of which (~ 69%) are located in the solar 
neighbourhood, within 3 kpc from the Sun. We estimate the age distribution of CSPN with average uncertainties of 
1-2 Gyr, and compare our results with the expected distribution based both on the observed mass distribution of white 
dwarfs and on the age distribution derived from available mass distributions of CSPN. 

Conclusions. We conclude most CSPN in the galactic disk have ages under 6 Gyr, and that the age distribution is peaked 
around 2-4 Gyr. 

Key words. ISM: planetary nebulae - Stars: AGB and post AGB - Stars: fundamental parameters 



1. Introduction 

Planetary nebulae (PN) are the offspring of intermediate 
mass stars with main sequence masses between 0.8 and 
8 Mq, approximately. As a consequence, their properties re- 
flect different physical conditions depending on the masses 
- and therefore ages - of their central stars (CSPN), which 
makes these objects extremely important in the study of the 
chemical evolution of the Galaxy (see for example Maciel 
et al. I2006|) . As an example, some recent theoretical models 
predict a time flattening of the observed radial abundance 
gradients in the galactic disk, while other models predict 
just the opposite behaviour. This can be analysed on the 
basis of abundance determinations in PN or in open clus- 
ter stars. In both cases, the results depend critically on the 
adopted ages of the objects considered. 

The determination of ages of CSPN is a complex prob- 
lem, and there is presently no single method that can be 
generally applied. In fact, most accurately determined ages 
refer to relatively young objects, for which methods such 
as lithium depletion, activity or cluster membership can be 
applied (see for example Hillebrand et al. (2009J) . Our group 
has pioneered in the treatment of this problem, and we 
have developed several methods to estimate the ages of the 
PN progenitor stars, based both on the observed nebular 
properties and in some properties of the stars themselves 
(cf. Maciel et al 120031 12005a| 12006} [2008]) . According to 
Soderblom (2009) , most age determination methods can be 
classified as (i) fundamental, (ii) model dependent, (iii) em- 



pirical, and (iv) statistical. The methods discussed in this 
paper belong to the last class. In principle, the traditional 
methods to derive the ages of galactic stars can be applied 
to CSPN, such as the use of theoretical isochrones (see for 
example Idiart et al. 2007), particularly for extragalactic 
nebulae. On the other hand, the physical properties of these 
objects are not as well known as in the case of normal stars, 
so that the derived isochrones are generally uncertain, lead- 
ing to the need of alternative methods. 

In this work, we discuss three methods developed so far 
by our group, namely the use of (i) an age-metallicity re- 
lation that depends on the galactocentric distance (Section 
2.1a), (ii) a simpler age-metallicity relation for the galac- 
tic disk (Section 2.1b), and (iii) the determination of ages 
from the central star masses obtained directly from the ob- 
served nitrogen abundance (Section 2.2). In Section 3, we 
estimate the expected age distribution of the CSPN based 
on (i) the observed mass distribution of white dwarf stars, 
and (ii) available mass distributions of CSPN, and compare 
the results with the distributions obtained by the methods 
mentioned above. In Section 4 a discussion is given, followed 
by some conclusions. 



2. Age determination and distribution of CSPN 
progenitors from nebular abundances 

2.1a. Method 1: The age-metallicity- galactocentric distance 
relation 
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Fig. 1. Age distribution of the central stars of planetary 
nebulae using an age-metallicity-galactocentric distance re- 
lation. The estimated age uncertainty is shown at the upper 
right corner. 
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Fig. 2. Age distribution of the central stars of planetary 
nebulae using the age-metallicity relation of Rocha-Pinto 
et al. (|2000p . The estimated age uncertainty is shown at 
the upper right corner. 



The first method to be considered was initially developed 
by Maciel ct al. (2003) , in the framework of an estimate 
of the time variation of the radial abundance gradients in 
the galactic disk. Using the oxygen abundance measured 
in the nebula, which is given in the usual form e(O) = 
log(0/H) + 12, the [O/H] abundance relative to the Sun 
is simply given by [O/H] = log(0/H) — log(O/H) , where 
we have adopted e(O) = 8.7 (see for example Asplund 
120031 and Asplund et al. 12003)) . The relation between the 
metallicity [Fe/H] and the oxygen abundance is given by 

[Fe/H] = 0.0317 + 1.4168 [O/H] , (1) 

which is valid for —1.0 < [Fe/H] < 0.5, as discussed by 
Maciel et al. (2003). Finally, the ages of the PN progenitor 
stars are given by an age-metallicity-galactocentric distance 
relation developed by Edvardsson et al. (1993), 

logt = 0.872 - 0.303 [Fe/H] - 0.038 R , (2) 

where t is in Gyr, and R is the galactocentric distance in 
kpc, so that some knowledge of the distance to the PN 
must be assumed. These relationships were applied to a 
sample of 234 well observed PN from Maciel et al. (|2003|) . 
located in the solar neighbourhood and in the galactic disk, 
for which the data are obtained with the highest accuracy. 
These objects have galactocentric distances in the range 
4 < i?(kpc) < 14, and most (69%) are located in the solar 
neighbourhood, with distances d < 3 kpc. The results are 
shown in Fig. 1. It can be seen that the distribution shows 
a prominent peak, located around 3-6 Gyr, similar to the 
age of the Sun, suggesting that most PN come from stars 
having masses close to one solar mass when in the Main 
Sequence. 

2.1b. Method 2: The age-metallicity relation for the galactic 
disk 

Rocha-Pinto et al. (|2000p derived an age-metallicity rela- 
tion for the galactic disk based on chromospheric ages and 
accurate metallicities. According to this relation (cf. Table 3 
of Rocha-Pinto et al. I200U|) . we can write 

[Fe/H] = 0.13969 - 0.08258 t + 0.00277 t 2 (3) 



where t is in Gyr. From this equation the stellar lifetimes 
can be determined once the metallicity is fixed. We can 
apply the same procedure as in the previous method (cf. 
Section 2.1a), and obtain [Fe/H] from the oxygen abun- 
dance e(O). The results for the same sample considered in 
Fig. 1 are shown in Fig. 2. The derived age distribution is 
flatter than in the previous method, but in both cases most 
stars have ages under 6 Gyr. 
2.1c. Uncertainty analysis 

Oxygen abundances are the best determined abundances 
in PN, with uncertainties typically under 0.2 dex, while 
other elements have generally higher uncertainties, in the 
0.2 to 0.3 dex range. In the case of the best determined 
O/H abundances, formal uncertainties under 0.10 dex are 
often estimated. The objects included in our sample are 
all disk planetary nebulae, avoiding the more distant ob- 
jects for which the abundances are poorly known. They re- 
sult from a very careful selection, in which the best spectra 
available were taken into account, and for which a compar- 
ison of abundances from different sources produces a very 
good agreement (see for example the individual abundance 
comparisons shown in Maciel et al. 2006"] ). All abundances 
considered are derived from collision excitation lines, which 
are considered as true representatives of the ionized gas in 
the nebulae (cf. Liu 12006). An additional uncertainty may 
be introduced by the ON conversion which may occur in 
the progenitor stars (cf. Stasihska 120081) . However, present 
results are not conclusive, and in any case they would only 
affect the progenitor stars near the upper mass bracket of 
the Main Sequence stars that produce planetary nebulae, 
which are a small fraction of the sample considered here. 

The solar oxygen abundance is accurate within 0.05 dex 
(see for example Asplund et al. I2006p . so that the uncer- 
tainty in the derived [O/H] abundances is essentially the 
same as in e(O). On the other hand, from the correlations 
presented in Maciel et al. (|2003p . iron and oxygen are clearly 
in lockstep in the galactic disk, so that we can safely adopt 
this hypothesis for our present sample. The average uncer- 
tainty in the [Fe/H] metallicity is dominated by the uncer- 
tainty in the [O/H] ratio, which is essentially the same as 
in the e(O) abundance, since the uncertainties in the slope 
and intercept of Eq. (1) are small, corresponding to 0.049 
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and 0.016, respectively. Therefore, an upper limit of about 
0.3 dex can be estimated for the uncertainty in [Fe/H], but 
a more realistic average would be about 0.2 dex, which cor- 
responds to an O/H uncertainty of roughly 0.2 dex. 

Concerning the ages as given by Eq.(2), a similar pro- 
cedure taking into account the [Fe/H] uncertainties leads 
to age uncertainties in the range 0.5 to 1 Gyr for objects 
in the solar neighbourhood, which are the majority in the 
sample considered here. 

The term involving the galactocentric distance provides 
a further uncertainty, but the range estimated above is 
not significantly changed. The abundances themselves arc 
distance-independent, but distance uncertainties may affect 
the galactic abundance gradients and as a secondary effect 
the position of the nebula projected onto the galactic plane. 
Independent confirmation of the reality of the gradients are 
presented by Peha et al. (|2005p . where a new method to 
determine the distances was developed, leading to a clear 
disk gradient that flattens out in the inner galaxy, a result 
confirmed by Pottasch and collaborators (cf. Gutenkunst 
et al. I2008[) . The main uncertainties in the gradients were 
discussed in detail by Maciel et al. (2005b), and the effect 
of the adopted distances was also analyzed by Maciel et al. 
(|2006|) and more recently by Maciel & Costa (|2009|) . In the 
latter work, we have taken into account four different PN 
distance scales, our own Basic Sample, and the scales by 
Cahn et al. (|TI)92") , Zhang (fT995|) . and Stanghellini et al. 
(2008). The main conclusion is that there is no appreciable 
change in the gradients when a sample of over 200 nebulae 
are considered with central stars in the age range of 2 to 10 
Gyr. It was also shown that, in average, the galactocentric 
distances obtained by these scales do not differ by more 
than about 1 kpc for objects in a ring centered at the so- 
lar position (Rq ~ 7.5 — 8.0 kpc) and extending to about 3 
kpc in either direction, which includes most nebulae in our 
sample. 

In agreement with these results, a typical uncertainty 
of cr(logi) ~ 0.10, where t is in Gyr, was estimated from 
the analysis by Edvardsson et al. ([1993), on which Eq. (2) 
is based. Therefore, the age distribution of Fig. 1 may be 
displaced by about 1 Gyr in either direction. Taking into 
account the observed scatter in the age-metallicity relation, 
this uncertainty is probably larger, up to about 0.2 dex (see 
a detailed discussion in Rocha-Pinto et al. 2000, 2006 ). As a 
consequence, absolute individual ages may have uncertain- 
ties higher than 1 Gyr, but it should be stressed that here 
we are interested in the age distribution, so that an actual 
displacement of the histogram of Fig. 1 by about 1 Gyr, as 
indicated by the horizontal bar at the upper right corner of 
the figure, is a realistic estimate. 

Considering now the uncertainties involved in the ap- 
proximation given by Eq.(3), we would like to stress that 
an age-metallicity relation is expected in the framework of 
a simple model of galactic chemical evolution (cf. Prantzos 
2008), and in fact several independent investigations have 
been able to derive some working relationships in the last 
two decades (cf. Rocha-Pinto et al. 120001 and references 
therein; Bensby et al. I2004)) . However, some problems are 
still under discussion, particularly the observed dispersion 
of this relation and its applicability to the thick and/or 
thin disks. The observed dispersion depends critically on 
the samples considered, and from the results by Evardsson 
et al. (|1993p and Feltzing et al. (|2001j) . a relatively large dis- 
persion of about 0.3 dex was estimated for the [Fe/H] ratio. 
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Our own results indicate that the actual dispersion may be 
considerably lower, about 0.2 dex (Rocha-Pinto et al.[2000, 
2006) , provided some important corrections are made, con- 
cerning the cosmic scatter, incompleteness of the samples 
and a careful analysis of stars that present contradictory 
age indicators, as discussed in detail by Rocha-Pinto et al. 
(2002). This value agrees very well with our estimate of 
the formal uncertainty, as discussed above. Therefore, the 
similarity of the uncertainties of the methods discussed in 
Sections 2.1a and 2.1b are reinforced, and the average age 
uncertainy of 1 Gyr is shown as a horizontal bar at the 
upper right corner of Fig. 2. These results are supported 
by a recent analysis by Pranztos (2()()Si. who presented 
a detailed discussion on the local age-metallicity relation 
based on the results of the Geneva-Copenhagen Survey (cf. 
Nordstrom et al. I2004|) . The resulting relationship shows 
a decrease in the metallicities with increasing ages, as ex- 
pected, up to about 6 Gyr, remaining essentially flat at later 
epochs. Taking into account the biases affecting this rela- 
tion by using simulated data as input for the age-metallicity 
relation, a monotonic relation was obtained, up to about 10 
Gyr ago, in better agreement with the views expressed in 
our previous work. 

2.2. Method 3: The age-N/O mass relation for CSPN 
2.2a. The age-N/O mass relation 

This method was also employed by Maciel et al. (|2003|) . and 
assumes a relationship between the central star mass mcs 
and the N/O abundance (for details see Maciel et al.[2003), 
which is expected from several theoretical and observational 
analyses. The adopted relation can be written as 

m CS = 0.7242 + 0.1742 log(N/0) (4) 

for -1.2 < log(N/0) < -0.26, and 

mcs = 0.825 + 0.936 log(N/0) + 1.439 [log(N/0)] 2 (5) 

for —0.26 < log(N/0) < 0.20. In these equations mcs is 
in solar masses. In order to obtain the stellar masses on 
the main sequence (tums), we adopted an initial mass- final 
mass relation of the form 

m C s = a + b m M s + c (m M s) 2 , (6) 

where a = 0.47778, b = 0.09028, anil r 0, to reproduce 
the same values in Maciel et al. (|2003jh namely mcs — 
0.55 Mq for m,MS — 0.8 M Q and mcs — 1-2 M for rriMS = 
8Mq, which reflect the known ranges of both CSPN and 
their progenitor stars on the main sequence. Concerning the 
mass-age relation, an approximation generally considered 
as accurate can be written as 

logi = d + e logrriMS + f (\ogm M s) 2 • (7) 

where t is the Main Sequence lifetime, or age, measured in 
Gyr, and the Main Sequence mass tums is given in solar 
masses. There are many discussions in the literature which 
suggest different values for the constants d, e, / (see for 
example Romano et al. I2005p . so that we will consider two 
cases here. In the first case, which we call Case A, we have 
d = 1.0, e = —2.0 and / = 0.0, which corresponds to a 
relation of the form t = C/m\ IS , where C — 10M|,Gyr. 
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Fig. 3. Age distribution of the central stars of planetary 
nebulae using a mass-N/0 abundance relation (Case A). 
The estimated age uncertainty is shown at the upper right 
corner. 

The second case, which we call Case B, adopts the well 
known mass-age relation by Bahcall & Piran (1983), which 
corresponds to d — 1.0, e = —3.6, and / = 1.0 in Eq. (7). 
Again we have t = 10 Gyr for tums = Taking into 

account Cases A and B has two advantages: first, it stresses 
that this important relation still involves some uncertain- 
ties; and second, the relation considered as Case B gives 
a steeper decrease of the stellar lifetimes at larger masses 
compared to Case A, so that we can easily interpret the be- 
haviour of the age distributions as a function of the mass- 
age relationship. 

Method 3 was applied to a sample of 122 PN for which 
all necessary data was available (cf. Maciel et al. [2003). 
Again we selected disk planetary nebulae for which the best 
abundance data were available so as to keep the uncertain- 
ties to a minimum. The results are shown in Figs. 3, and 4 
for cases A and B, respectively. These results are similar to 
the age-metalllicity-radius method, in the sense that most 
objects have ages lower than about 10 Gyr, and there is a 
sharp maximum in the probability distribution. However, 
its location depends on the estimated lifetimes as a func- 
tion of the Main Sequence mass, being around 3-7 Gyr for 
Case A and 1-4 Gyr for Case B. From case A to case B the 
lifetimes of the more massive stars are decreased, so that 
the probability of finding stars at larger lifetimes decreases, 
and the whole peak moves to the left, as shown by Figs. 3 
and 4. 

2.2b. Uncertainty analysis 

An estimate of the uncertainties involved in Method 3 sug- 
gests that they are similar or somewhat larger than in 
Methods 1 and 2. The basic relation is given by Eqs.(4) and 
(5), which relate the central star masses and the N/O abun- 
dances. Although such a relation is expected from theoret- 
ical models (see for example Renzini & Voli 1981, Marigo 
IMKH |2"00"TI Perinotto et al. IMHj) , there is presently no 
clearcut functional dependence between these quantities, 
which led Cazetta & Maciel ({2000)) and Maciel (|2000|) to 
propose a calibration based on the best parameters avail- 
able. The N/O ratio is relatively well determined, with un- 
certainties similar to the O/H ratio or even better, namely, 
under 0.2 dex, as already discussed. Adopting this value 
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Fig. 4. The same as Fig. 3 for Case B. 



and also taking into account the uncertainties in the co- 
efficients of Eqs. (4) and (5), an upper limit of about 0.2 
Mq would be obtained for the mass uncertainty. However, 
a more realistic estimate would be lower than this value, 
as suggested by the recent application of N/O masses to a 
well studied sample of CSPN by Maciel et al. (|2lK?g]) . 

In that work, the derived masses were successfully com- 
pared with the known mass distributions of CSPN and 
white dwarf stars in order to explain the relationship be- 
tween the modified momentum of the CSPN winds and 
the stellar luminosity. Some additional support to the N/O 
masses comes from the reanalysis of Tinkler & Lamers 
(2002| of the central star masses determined by Kudritzki 
et al. (|1997p . These results were based on an homogeneous 
set of parameters obtained from Zanstra temperatures, dy- 
namical ages and evolutionary tracks. The average masses 
obtained by Tinkler & Lamers (|2002p for the same sample 
studied by Maciel et al. (20080 and the N/O masses are in 
very good agreement. In view of these considerations, an 
average uncertainty closer to 0.1 Mq would be more appro- 
priate to Eqs. (4) and (5). 

Considering now the initial mass-final mass relation as 
given by Eq.(6), there are many determinations of this re- 
lation in the literature (Vassiliadis & Wood [T5551 Marigo 
mm mm Marigo & Girardi IMJTl Meng et al. MM, 
which generally agree with each other within 0.1 to 0.2 Mq, 
so that the total uncertainties in the stellar masses are prob- 
ably not much affected by this relation. As a consequence, 
the average uncertainties in the indivdual CSPN ages cor- 
respond to a mass uncertainty of at most 0.2 Mq. In fact, 
the main contribution to the age uncertainty is due to the 
age estimates as given by Eq.(7), and at this point the best 
procedure to overcome such difficulty is to adopt some kind 
of parametrization, by considering cases A and B. Applying 
the mass uncertainty mentioned above, a formal age uncer- 
tainty would be about 3.1 and 4.7 Gyr for cases A and B, 
respectively, considering a typical star of one solar mass 
at the Main Sequence. These value can be considered as 
upper limits, as they were obtained using a 0.2 Mq mass 
uncertainty. Adopting the more reasonable value of about 
0.1 Mq, the average uncertainties would be 1.7 Gyr and 2.9 
Gyr for Cases A and B, as shown by the error bars in Figs. 3 
and 4. Again, it should be stressed that we are interested 
in the age distribution, rather than in individual ages, so 
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that the general effect of the uncertainties is expected to 
be smaller than indicated by the formal uncertainties. 

2.2c. Binary CSPN 

A final comment is appropriate on Method 3, concerning 
the possibility that a significant fraction of the PN may 
have originated from close binary stars. Current results 
are controversial, in the sense that a very small number 
of binary PN are known, as discussed by De Marco (2009), 
while some investigations suggest that binaries constitute a 
much larger proportion of the galactic nebulae (see for ex- 
ample Mi szalsk i et al. I2"00"9l De Marco l2"0"0l)I [MM Moe & 
De Marco [21305)) . In view of these contradictory aspects, it is 
difficult to evaluate the effect of binarity on the present re- 
sults. It would be expected that the age distribution would 
be displaced towards higher ages by only a small amount, 
as suggested by the good agreement of the results of our 
methods with the actual mass distribution of CSPN and 
white dwarf stars, as we will see in the next section. 



3. The expected age distribution of the central 
stars of planetary nebulae 

3.1. The white dwarf mass distribution 

The expected age distribution of the CSPN can be esti- 
mated by analyzing the much better known mass distribu- 
tion of the white dwarf stars. Since the average mass loss 
rates during the PN phase amount to about dM/dt ~ 10~ 8 
to 10~ 6 M Q /yr and the whole PN phase duration is about 
At ~ 1 to 2 x 10 4 yr, the total mass lost during this phase 
is Am ~ (dM/dt) At ~ 10~ 4 to 2 x 1O _2 M , which is 
much smaller than the CSPN masses. Therefore, as a first 
approximation, the mass distribution of the CSPN must be 
similar to that of the white dwarfs, except for the very low 
mass stars with m < 0.55 Mq. Such stars are not expected 
from theoretical models, since main sequence stars lead- 
ing to white dwarfs with masses lower than about 0.55 Mq 
probably go directly to the white dwarf phase. 

Recent work on the mass distribution of white dwarfs 
by Madej et al. (|2"DM|) and Kepler et al. (|2"0"0"7) lead to a 
distribution which is strongly peaked at about 0.56 M©, as 
shown in Fig. 5 (cf. Madej et al. I2004)) . This investigation 
was based on a large sample of about 1200 white dwarfs, 
and shows very few objects (about 5%) with masses larger 
than about 0.8 Mq. 

The white dwarf mass distribution can be well fitted by 
a Gaussian probability density distribution defined by 



. , N(m) 1 
P( m ) = -^—L - — 

y ' N t Am 



(8) 



where N(m) is the number of stars with mass m, Nt is the 
total number of stars, and Am = 0.025 Mq is the size of 
the adopted mass bins, as shown in Fig. 5, where the curve 
plotted is 30 x P(m). As can be seen, the gaussian fit to 
the data is good, although at very small and very large 
masses the predictions are somewhat lower than observed. 
Such probability is normalized, that is, J P(m) dm = 1. 
For white dwarfs P(m) can then be written as 




m(wd) (M ) 

Fig. 5. Mass distribution of white dwarfs (Madej et al. 
I2004P and Gaussian fit. 




m(MS) (M ) 

Fig. 6. Mass distribution of CSPN progenitor stars at the 
Main Sequence. 



where A = 0.90349, a = 0.08798 and \l = 0.55420. 
Assuming that the CSPN have approximately the same 
mass distribution than the white dwarfs, we have 
Pcs( m ) dm = Pwr){m) dm. As a first approximation, we 
can assume a linear relation between the CSPN mass and 
the Main Sequence (MS) mass, as given by Eq. (6) with 
c = 0, so that we have then a = 0.47778 and b = 0.09028, 
as before. The main sequence masses can be written as 



mMS = t (mcs - a) 



(10) 



P{m) 



A 



/2tt 



(9) 



and dmcs — b drriMSi so that the probability distribution 
for MS stars can be written as PMs( m ) = b Pcs(m). The 
derived probability distribution of MS stars is shown in 
Fig. 6. It can be seen that it peaks around one solar mass, 
and also includes some very low (even negative) masses, 
which of course are unrealistic, as they correspond to CSPN 
masses lower than mcs — a — 0.48 Mq, which are not ob- 
served. However, as we will see later on, this essentially 
increases the probability of very high lifetimes, an effect 
that can be interpreted with some modifications in the cal- 
culation of the probability PMs( m )- 

Assuming that the star formation rate in the Galaxy has 
remained approximately constant along the galactic life- 
time, the age distribution of the CSPN can be estimated 
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Fig. 7. Age distribution of the progenitors of the central 
stars of planetary nebulae. Solid line: case A, dashed line: 
Case B. 



from the mass distribution of their progenitor stars. In or- 
der to derive the age distribution of the observed CSPN, we 
will adopt cases A and B considered previously. Considering 
that P(t) dt = PMs(m) dm, we have for case A: 



P(t) = ^ Pusim) 



and for case B: 



P(t) = 



,4.6 



2 log to- 3.6| lfji+Oog™) 2 



PMs{m) 



(11) 



(12) 



where we have dropped the subscript MS from the mass. 
The obtained age distributions are shown in Fig. 7 for Case 
A (solid curve) and B (dashed curve), respectively. Case A 
gives larger lifetimes for masses greater than one solar mass, 
which pushes the peak of the probability distribution to the 
right (2-3 Gyr), while for case B these lifetimes are shorter, 
and the peak moves to the left (~ 1 Gyr). Excluding the 
MS stars that do not lead to the formation of CSPN, the 
main effect on Fig. 7 is a sharp decrease in the probability 
for ages greater than about 12 Gyr, leaving the peak re- 
gion essentially unaffected. The main conclusion that can 
be drawn is that a peaked distribution can be expected, but 
the precise location of the peak depends on the adopted ini- 
tial mass-final mass relation and especially the calculation 
of the lifetimes. 

For the sake of generality, we have also considered a 
quadratic relation as in Eq. (6) with c ^ instead of a 
linear equation. In this case, we have adopted the con- 
stants a, b, and c from the recent work of Meng et al. 
(|2008|) . taking the metal abundance Z = 0.02. The cen- 
tral star mass is given by the minimum of mcsipi, bi, c i) 
and TOcs(a2, &2, C2), where 01 = 0.5716, b\ — —0.04633, 
ci = 0.02878, a 2 = 1.1533, b 2 = -0.2422, and c 2 = 0.04091. 
In this case, the probability Fms(w) is given by PMs( m ) — 
(6+2 c Toms) Pcs( m )- As it turned out, the results of Fig. 7 
are not particularly sensitive to this assumption. 

3.2. The observed mass distribution of the central stars of 
planetary nebulae 

The mass distributions of CSPN and white dwarfs have 
also been previously considered by Stasihska et al. (|1997|) 



and Napiwotzki (2006). More recently, Gesicki & Zijlstra 
(2007) analyzed these distributions based on a dynamical 
method which allows mass determinations within 0.02 M 
(cf. Gesicki et al. [2006). The CSPN masses were obtained 
for a sample of 101 objects from a relation between the 
temperatures of the central stars and the dynamical age 
of the surrounding nebulae. The ages were derived from 
a combination of recent spectra, line ratios and nebular 
sizes, using a photoionization model to obtain the central 
star temperature. Theoretical tracks by Blocker (19951) were 
used to derive the stellar masses. It results that both the 
CSPN and white dwarf distributions peak around 0.6 Mq as 
in Madej et al. (|2004p and Stasihska et al. (|1997p . although 
the white dwarf distribution shows a broader mass range. 
The CSPN distribution shows essentially no objects with 
masses higher than 0.7 Mq, while in the case of the DA 
white dwarfs, which are presumably the offspring of H-rich 
CSPN, a few objects are observed with masses larger than 
0.8 Mq. Although there may be some differences between 
the recently obtained white dwarf mass distributions, as 
discussed by Gesicki & Zijlstra (2007)), they all agree in the 
sense that any sizable sample of CSPN is expected to have 
a larger number of objects with masses close to 0.6 Mq. 
These results are in good agreement with our own N/O 
masses, as discussed by Maciel et al. (|2008p . 

Instead of adopting the white dwarf mass distribution, 
we may then use directly the mass distribution of CSPN as 
recently derived by Gesicki and Zijlstra (20071). The mass 
distribution is shown in Fig. 8, and can also be approxi- 
mated as a Gaussian probability density distribution simi- 
lar to eq. (9), where we find A = 0.89627, o = 0.01431, and 
fi = 0.6087. The corresponding probability distribution for 
CSPN is also shown in Fig. 8. Adopting again the same 
hypotheses as in Section 3.1, the CSPN age distributions 
are shown in Fig. 9 for Case A (solid curve) and B (dashed 
curve), respectively. In this case, the peaks are located at 
approximately 4-6 Gyr (Case A) and 2-4 Gyr (Case B). 

4. Discussion 

Some indication of the age distribution of CSPN may be 
obtained from the PN classification originally proposed by 
Peirnbert (|1976[ H"990p . According to this classification, the 
following approximate ages are attributed to the different 
PN types: Type I (1 Gyr), Type II (3 Gyr), Type III (6 
Gyr), and Type IV (10 Gyr) (see for example Stasihska 
2004, table 7). Type IV are halo objects, which are rarer and 
presumably much older than the remaining types, which are 
located in the disk and bulge of the Galaxy. Since nebulae 
of Types I-III constitute the vast majority of the known 
galactic PN, in this work we will not take into account the 
older, Type IV nebulae. 

Fig. 7 shows the CSPN age distributions for Cases A 
and B using the mass distribution of white dwarfs (Madej 
et al. I2004p . and Fig. 9 gives the corresponding distribu- 
tions from the CSPN (Gesicki and Zijlstra 2007) mass dis- 
tribution. The merged age distributions suggest a maxi- 
mum around 2-6 Gyr for Case A and 1-4 Gyr for case B. 
Again, the main conclusion that can be drawn is that a 
peaked distribution can be expected, but the location of 
the peak depends on the adopted assumptions. In prin- 
ciple, we would expect Case B to be more realistic than 
Case A (see the discussion on stellar lifetimes by Romano 
et al. I2005p . so that the results of Section 3 would suggest 
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Fig. 8. Mass distribution of the central stars of planetary 
nebulae from Gesicki & Zijlstra (|2007|) and Gaussian fit. 

a preferable range of 1-4 Gyr for the peak of the distribu- 
tion. On the other hand, the results based on the empirical 
CSPN mass distribution are probably more accurate than 
those inferred from the white dwarf mass distribution, for 
the reasons mentioned in Section 3. Therefore, in view of 
these considerations, and taking into account the uncer- 
tainty analyses of Sect. 2.1c. and 2.2b., we suggest that the 
peak of the age distribution is probably located around 2-4 
Gyr, as shown by the dashed curve of Fig. 9. We may then 
compare this distribution with the results of Sect. 2. Taking 
into account the uncertainties involved in the age determi- 
nations, which are estimated to be in the range 1.0-2.0 Gyr, 
approximately, as discussed in Sect. 2, it can be concluded 
that all methods considered produce results reasonably in 
agreement with the expected age distribution. According 
to the discussion in Sect. 2, Methods 1 and 2 have similar 
results, as indicated by Figs. 1 and 2, but the age distribu- 
tion of Method 2 is less strongly peaked than either of the 
distributions obtained in Sect. 3 (cf. Figs. 7 and 9). Since 
the uncertainties are similar, Method 1 is probably more 
accurate than Method 2. 

At face value, the age distributions of Method 3 shown 
in Figs. 3 and 4 are very similar to the results of Fig. 9, 
especially if Case B is considered, as seems more appro- 
priate. In other words, the best results are probably those 
by Method 3, Case B (Fig. 4), which show a very good 
agreement with the results from the empirical mass distri- 
bution of CSPN, as shown by the dashed curve of Fig. 9. 
Although the estimated uncertainty of Method 3 is larger, 
it is reassuring that the average distribution shows such a 
remarkable agreement with the expected age distribution 
of CSPN. Therefore, our results obtained from completeley 
independent methods and sources of data are reasonably in 
agreement with each other, so that we can have an estimate 
of the age distribution of CSPN. Naturally, the details of 
the individual age determinations still need to be worked 
out. 
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